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In light of the experimental evidence for the existence of a superflnidity region in a Bose-condensed 
sodium gas a theoretical model is put forward. It will be shown that the predictions of the present 
work do match with the extant measurement readouts. As a byproduct we also calculate the speed 
of sound and compare it against the current experimental results. 
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The first experimental evidence concerning the phe¬ 
nomenon of superfluidity can be tracked down to the 
work of Kamerling Onnes in the year of 1911 with helium 
when it was found that if cooled below 2.2 °K He did not 
contract but rather expand [1]. Since then the amount 
of theoretical and experimental work has been able to 
provide a coherent picture to the subjacent Physics [2- 
4]. Another phenomenon, which also emerged in the last 
century, is the Bose-Einstein Condensation (BEC) the 
one, as in the case of superfluidity, is inherently related 
to the presence of very low temperatures. The need for 
low temperatures for the appearance of these two effects 
leads us to the obvious question concerning a possible 
relation between them. It was Fritz London [5] who put 
forward the idea of a connection between these two effects 
stating that the transition from He I (the high temper¬ 
ature phase of liquid helium) and Hell (the low tem¬ 
perature phase) should be considered an example of a 
BEC. Even more, he suggested that in He II a macro¬ 
scopic quantum current of matter could be present, i.e., 
he introduced the idea that BEC and superfluidity could 
appear, simultaneously, in a system. 

It has to be stressed that this understanding cannot 
be considered a closed issue [6], and the answer to this 
aforementioned interrogant is that, though there is a 
close relationship, it is not a unique one. Indeed, we 
may state that BEC is neither necessary nor sufficient 
for the presence of superfluidity. For instance, and ideal 
Bose-Einstein condensate shows no superfluidity, and, as 
a counterpart a two-dimensional superfluid cannot con¬ 
densate [7]. 

Landau (within the two-fluid model proposed by Tisza 
[8]) introduced the concept of elementary excitations [9] 
as a fundamental element in the description of the be¬ 
havior of He II. Landau argued [10] that the normal 
fluid (the non-superfluid component) could be regarded 
as a dilute gas whose components are weakly-interacting 
elementary excitations which move in a background de¬ 
fined by the superfluid component. Accordingly, the phe¬ 
nomenon of superfluidity appears if the velocity of the 
corresponding flow lies below a certain threshold value 
given by 
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( 1 ) 


Here e{p) denotes the energy of an elementary excita¬ 
tion and p its corresponding momentum. If the velocity 
is larger, then the microscopic rugosities of the walls of 
the container will scatter the particles of the fluid entail¬ 
ing the loss of kinetic energy of the fluid, i.e., viscosity 
appears. 

In the experimental realm the quest for this critical 
velocity has been carried out in a sodium-BEC, and the 
results show a possible velocity threshold located around 
the value of 1.6mm/s [II]. The use of sodium-condensed 
gas in the experimental context is related not only to the 
aforementioned case but also to the excitation of phonons 
by light scattering [12] and the propagation of sound [13, 

14]- 

The main purpose of the present work is to obtain a 
prediction for the critical velocity for a BEC. The de¬ 
duced model will be compared against the reported mea¬ 
surement readouts [11], and, in addition, the speed of 
sound for a sodium-condensed gas will be found and com¬ 
pared with the current experimental results [13, 14]. 

From a fundamental point of view our mathematical 
model can be defined by an A^-particle Hamiltonian the 
one in the formalism of second quantization is [7] 


H= dr 
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In this Hamiltonian and tjj{r,t) represent 

bosonic creation and annihilation operators, respectively. 
It is restricted to low energies and momenta and implies, 
as a consequence of the aforementioned conditions, that 
the interaction among the particles is, as usual, codified 
by the scattering length parameter a, i.e., Uo = . 

The trapping potential V (r), for our case, corresponds to 
an isotropic harmonic oscillator whose frequency reads to. 
In addition, there are N particles in the gas, each of them 
with mass m, the volume occupied by the system is V. 
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Our mathematical assumptions are: 

(i) Only the ground and the first excited states are 
populated. This condition can be justified recalling that 
for a bosonic system, with chemical potential /i and en¬ 
ergy levels of single-particle e, the occupation number in 
thermal equilibrium is given by [15] (j3 = 1 /{kT)) 


Clearly, 

N(o)= j {'ip(o){r)Yd^r, (9) 


< >= eC-i/^-1 - 

It is readily seen that we deal with a monotonic de¬ 
creasing function of e, and this feature justifies the 
present assumption. 

(ii) The mathematical description of the two occupied 
states will be done resorting to the Hartree approxima¬ 
tion, in which the ground state of the interacting system 
is deduced by a Ginzburg-Pitaevski-Gross energy func¬ 
tional [16], and it entails that the ground state wave- 
function corresponds to the case of a harmonic oscillator 
situation but the frequency is modified due to the fact 
that the system has a non-vanishing scattering length 
[17], such that the fundamental length parameter reads. 

Here I is the radius related to the trap given by the 
isotropic harmonic oscillator, namely. 



( 10 ) 


The wavefunction of the first excited state will be con¬ 
sidered as the first excited state of an isotropic oscillator 
related to frequency given by (6) and, due to our sym¬ 
metry, it has three possibilities, all with the same math¬ 
ematical structure, namely. 
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Here = x, = y, and x^^'> = z. 

Of course, (11) must be related to the total number of 

particles in excited states (fV(e) = |^), a condition 
that becomes [15] 


Nie) = 




2=1 


d-^r. 


( 12 ) 



Clearly this condition implies an effective frequency 


Having stated our assumptions we proceed to compute 
the speed of sound and the critical velocity. The energy 
of the ground state in its three possibilities, i.e., kinetic 
energy, due to the trap, and interaction are [15, 17] 
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Usucilly tli6 6xp6riiii6ntcLl conditions cntnil ^ I [11] 
and, in consequence, u) < uj. 

In other words, the order parameter related to the par¬ 
ticles in the ground state is provided by 


J V(f)'ip(^o){r)d^r = ^TOW^i?^iV(o). (14) 
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In this last expression denotes the number of par¬ 
ticles in the lowest energy state. The presence of a non¬ 
vanishing scattering length entails that in the ground 
state not all the particles can have zero-momentum, the 
reason for this lies in the fact that the two-body interac¬ 
tion mixes in components with atoms in other states [15] 
and 


N(o) 


= N 
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The energy of the ground state (here denoted by if(o)), 
no elementary excitations are present, is the sum of 
the last three expressions. The corresponding pressure 

{P{o) = -dE(Q)/dV) is 
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If z;(o) denotes the speed of sound related to the last 
expression = —{V‘^/mN)dP(^Q'^/dV) we are led to 
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For the case of sodium [13, 14] the physical param¬ 
eters are: m = 36.8 x 10~^'^kg, a = 2.75 x 
N = 5 X 10®. Considering one of the reported peak 
densities of the condensate, namely, 1 x 10^^cm“®, we 
obtain, in the roughest approximation from our calcu¬ 
lations, 6.59 mm/s. Clearly, our prediction is in good 
agreement with the experiment [13]. 

We now proceed to compute the lowest energy and mo¬ 
mentum of the elementary excitations, physical parame¬ 
ters required for the deduction of the critical velocity [3]. 
The deduction of the energy of an elementary excitation 
and of its corresponding momentum requires the knowl¬ 
edge of the energy of a single-particle in the first excited 
state [7]. Our assumptions entail that the thermal cloud 
contains particles subject to an isotropic harmonic oscil¬ 
lator whose frequency is (6) therefore the energy of an 
excited particle is given by this assumption and easily 
calculated as a function of the effective frequency of our 
variational procedure 


the occupation number of the particles in the first ex¬ 
cited state fulfill the condition A^(g) ~ 10^ > 1, and, 
in consequence, < >= 1. Our assumptions im¬ 

ply < hgjp >= 0, Vi > 1. Indeed, we have consid¬ 
ered that the thermal cloud is comprised by particles 
which occupy only the first excited state, in other words, 
< ngj.j >= 0, Vi > 1. Introducing this condition into 
(21) leads us to the aforementioned result for the occu¬ 
pation number of the elementary excitations. 

Casting (19) in terms of the effective volume V = 
47ri?®/3, and using (4), (5), and (6), we have that the 
energy of our elementary excitation is 


e 
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We must now find the momentum of this elementary 
excitation. Elementary excitations, which define the nor¬ 
mal component of the fluid fluid, can be regarded as 
a bosonic gas whose components are weakly-interacting 
and moving in a region in which a constant potential ex¬ 
ists, and this potential is defined by a mean field approach 
[7]. According to this interpretation we may rewrite (22) 
in the same form as in the case in which our BEC is a 
homogeneous one [7]. In other words, we cast our last 
expression in the following form 


e 


2m\l ^ V 


(23) 


(18) 

According to Bogoliubov [7, 18] the energy of an ele¬ 
mentary excitation, here denoted by e, is a function of 
the energy of the excited particles of the BEC, namely. 


e = 



2NU(o) - 


(19) 


The energy of all the elementary excitations turns out 
to be [7, 18] 


E = ( 20 ) 

Here < hg > denotes the occupation number of the ele¬ 
mentary excitations with energy e. The relation between 
the occupation numbers of particles and elementary ex¬ 
citations is [7] 


< fie >= 


< He > 

1+ < He > 


( 21 ) 


At this point, for the sake of simplicity, we resort to 
the experimental values related to the detection of a crit¬ 
ical velocity in a sodium condensed gas [11] in which 


Clearly, (23) allows us to deduce the wavenumber re¬ 
lated to our elementary excitation and, in consequence, 
its momentum. Indeed, we have for these two physical 
variables, respectively, that 
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Resorting to Landau criterion (1) we obtain that the 
critical velocity is given by 

1 h 125 /47r\ 2/3 207riVa 

U(gHt) = +-^7T7r- (26) 

The experimental parameters [11] are a critical speed 
of = 1.6 mmjs. In addition, the number of parti¬ 

cles in this experiment has a minimum of TV = 3 x 10® 
and a maximum of = 12 x 10®, and for the evaluation 
of our expression we will take the arithmetic average, i.e., 
N = 7.5 X 10®. The effective volume is that of an ellipsoid 
whose axes are /i = 45 x 10“®m and li = 150 x 10“®m 
such that V = 
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Introducing these values into (26) entails 


y(m) _ 2 ^ gg mmjs. 

(27) 

The reported critical speed is [11] 


= 1.6 mmjs 

(28) 

The ensuing error is less that 18 percent 


1^,0 = 0.179. 

(29) 


Since the number of particles in the corresponding ex¬ 
periment varies from N = 3xlO®to-/V = 12x10® 


[11] the associated values for the critical speed go from 
1.24 mmjs to 2.48 mmjs. If iV = 5 x 10®, then 

yim) _ y{e) ^ 


In conclusion, we have put forward a theoretical model 
for the deduction of the critical velocity in a sodium- 
condensed gas. This threshold speed has been computed 
and compared against the extant experimental results, 
having a good agreement between them. A more precise 
evaluation of the present idea requires a better knowledge 
of the value of N employed in the experiment. Previous 
works offer larger critical velocities which have a bigger 
error than the one here deduced [11, 19], when compared 
to the experimental result. 
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